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Abstract 

We describe the spin distribution in the vicinity of a non-magnetic impurity in a two-dimensional 
antiferromagnet undergoing a transition from a magnetically ordered Neel state to a paramagnet 
with a spin gap. The quantum critical ground state in a finite system has total spin S = 1/2 (if the 
system without the impurity had an even number of 5 = 1/2 spins), and recent numerical studies 
in a double layer antiferromagnet (K. H. Hdglund et al, Phys. Rev. Lett. 98, 087203 (2007)) have 
shown that the spin has a universal spatial form delocalized across the entire sample. We present 
the field theory describing the uniform and staggered magnetizations in this spin texture for two 
classes of antiferromagnets: (i) the transition from a Neel state to a paramagnet with local spin 
singlets, in models with an even number of 5 = 1/2 spins per unit cell, which are described by a 
0(3) Landau-Ginzburg- Wilson field theory; and (ii) the transition from a Neel state to a valence 
bond solid, in antiferromagnets with a single 5 = 1/2 spin per unit cell, which are described by a 
"deconfined" field theory of spinous. 
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I. INTRODUCTION 



There have been many experimental studies of non-magnetic Zn impurities substituting 
for the spin 5 = 1/2 Cu ions in spin-gap and superconducting compoundsi'^i'^'^'^'^. These 
have stimulated many theoretical studies of the spin dynamics in the vicinity of a vacancy 
{i.e. a site with no spin) in S" = 1/2 square lattice antiferromagnets2.i^i2ii2iiiii2iii2iiiiii^ii^iiZii^ii£. 

An important feature of the impurity-response escaped^ theoretical attention until 
recently^. Consider the regime where the bulk antiferromagnet preserves global ro- 
tational symmetry and has a S = ground state. Such states can be reached by 
deforming the nearest-neighbor antiferromagnet into a coupled-ladder or coupled-dimer 
antiferromagnet^"^, in a double-layer antiferromagnet^, or by adding additional ring- 
exchange interactions while preserving full square lattice symmetry2^. Now remove a single 
S" = 1/2 spin in a system with an even number of spins, leaving an antiferromagnet with a 
vacancy and an odd number of S = 1/2 spins. We expect this antiferromagnet to have a 
doubly-degenerate ground state with total spin S = 1/2. Without loss of generahty, we can 
examine the ground state with spin-projection 5"^ = 1/2. In such a state, even though there 
is no broken symmetry and no applied magnetic field (the Hamiltonian has full SU(2) spin 
symmetry), the expectation values of the spin projection on the site i, (Szi), is non-zero on 
all i for any finite system of size L. The question of interest in this paper is the following: 
What is the spatial form of {Szi) ? It is possible that the = 1/2 magnetization is pushed 
out to the boundaries of the system, far from the impurity: in this case, it will not be rel- 
evant to the impurity properties in the limit L ^ oo. However, we will find this is not the 
case for the antiferromagnets examined in this paper. For the spin-gap antiferromagnets we 
consider, the S" = 1/2 magnetization is bound to the impurity over a length scale inversely 
proportional to the spin gap. At the quantum critical points separating the spingap states 



from the Neel state, which define 'algebraic spin liquids', we will find, as in Ref. [181, that the 
impurity magnetization is delocalized over the entire system, forming a spin texture with a 
universal spatial form determined only by the system size L. 



First, let us consider the models which have been numerically studied in Ref. Il8|. These 
are antiferromagnets which have an even number of = 1/2 spins per unit cell (such as the 
coupled-dimer— 1^ or double layer— models), which exhibit a transition between a Neel state 
and a simple spin gap state; the latter state is adiabatically connected to a state in which 
the spins in each unit cell are separately locked into singlets, with negligible resonance 
between unit cells. This is a 'conventional' transition, described by a Landau-Ginzburg- 
Wilson (LGW) theory. A convenient description of both phases and the quantum phase 
transition is provided by the 0(3) non-linear sigma model, expressed in terms of a unit 
vector field n(x, r) representing the local orientation of the Neel order parameter. Here 
X is the two-dimensional spatial position, r is imaginary time, and = 1 everywhere in 
spacetime. The bulk action in the absence of the impurity is the 0(3) non-linear sigma 



^ Section II. B. 3 of Ref. 115! contains results which can be used to extract the spin textures in zero field. 
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model 

5b =^/rf^/ dMd,n)\ (1.1) 

where g is the couphng constant which tunes the antiferromagnet from the Neel state {g < gc) 
to the spin gap state {g > gc), fi is a 3-dimensional spacetime index and a spin- wave velocity 
has been set to unity. In this formulation, the influence of the impurity is represented 
universally by the following Berry phase term alone^^ (provided the antiferromagnet is not 
too far from the critical point) 

^r^p = ^SJ dr A[n{0, r)] ■ (1.2) 

for a spin S = 1/2 antiferromagnet, where A is the Dirac monopole function in spin space 
with Vn X A = n. Note that does not include any coupling constants, and it depends 
upon the value of n only at :r = 0, which is the position of the impurity. 

Now we need to describe the 5 = 1/2 ground state of S^ + S^^ for g > gc- First, we need 
a proper discussion of the rotationally invariant 5 = ground state without the impurity. 
While it may be possible to do this within the context of a small g expansion of the 0(3) 
non-linear sigma model, the procedure is quite cumbersome and delicate, requiring a global 
average over all possible locally ordered states. We shall instead follow a simpler procedure 
which is described in more detail in Section [III we use an alternative soft-spin, LGW for- 
mulation of iS" in terms of a vector order parameter, (p, whose length is unconstrained. The 
= saddle point then is an appropriate starting point for describing the physics of the 
S = ground state of the bulk theory and its excitations. Next, we include the impurity 
term described by S^^, and also apply an infinitesimal magnetic field in the z direction. As 
we will show in Section |Tll the Berry phase effectively localizes the order parameter at the 
impurity site, n(x = 0, r), to a specific orientation on the unit sphere; in particular, for the 
Sz = 1/2 state chosen by the apphed field, we may perform an expansion about a saddle 
point with n(x = 0, r) = (1, 0, 0). This expansion quantizes, at each order, the total spin at 



5*2 = 1/2: this was established in Section 11. C. 2 of Ref. Il0| for g < g^ and the same result 
also applies here for g > gc- The infinitesimal magnetic field is set to zero at the end, but 
the spin density of the = 1/2 state remains non-zero in this limit. 

The results in Section |TT] provide an explicit analytic realization for the scaling forms 



presented in Ref. Il8| for the spin texture near the impurity. For the magnetization density, 
Q, which is the conserved Noether "charge" density associated with the 0(3) symmetry of 
the antiferromagnet, we have ai g = gc and zero temperature (T) and in the Sz = 1/2 state: 

{QM) = ^^Q (t) (1-3) 



L2 \L 
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where is a universal function obeying the quantized total spin condition 



S. 



[lA) 



Similarly, the staggered magnetization associate with the Neel order parameter obeys the 
scaling form 



1 



L(i+^)/2 



X 



:i.5) 



a.t g = Qc, where $„(r) is another universal function, but its overall scale is non-universal. 
The exponent rj is the anomalous dimension of n at g = Qc in the absence of the impurity. 

Let us now turn to the more interesting and much more subtle case of a "deconfined" 
critical point^i^. Here we are considering antiferromagnets with a single S = 1/2 spin per 
unit cell, and so there is no simple spin-gap state with local singlets. For the models studied 
in Refs. l23lJ24j]25l . the spin gap state has singlet valence bonds which crystallize into a regular 



arrangement, breaking the space group symmetry of the square lattice, while preserving spin 
rotation invariance. Such a state is a valence bond solid (VBS), and we are now interested in 
the impurity response across the Neel- VBS transition. As argued in Refs. l24lJ25l the vicinity 
of this quantum critical point is described by the CP^ field theory which is expressed in terms 
of a bosonic 'spinon' represented by a complex spinor field Za{x, r), where a =t, |, and the 
constraint l^al"^ = 1 is obeyed everywhere in spacetime. The Neel order parameter, n is 
related to by 

n = zlaapzp, (1.6) 



where a are the Pauli matrices. Also, in our analysis, we find it useful to generalize to 
the CP^^^ model with SU(A^) symmetry, where a = 1 . . . A^, and then the Pauli matrices 
are replaced by the generators of SU(A^). The action of the CP^~^ model also involves a 
non-compact U(l) gauge field A^, and is given by 



dr d X 



\{d^-tA^)z^\^ + ^' 



:i.7) 



This theory describes a Neel-ordered phase for g < g^, and a spin-gap state with VBS 
order for g > gc (additional Berry phase terms are needed to obtain the four-fold square- 
lattice symmetry of the VBS order-^). It is crucial to note that, unlike the situation in 1-1-1 
dimensions^i^^, the models iS" (in Eq. fll.ll)) a nd are not equivalent to each other in 2-1-1 
dimensions. This was established in Ref. |29!, and is a consequence of the proliferation of 
'hedgehog' or 'monopole' defects at the critical point of 5"; such defects are absent in the 



theory. 
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that 



Now let us add an impurity to the field theory in Eq. (11.71) . It was argued in Ref. 
the impurity is now represented by a source term for a static charge Q = 2S at x = 0. Thus 



C2 



iQ / drArix = 0, r) 



(1.8) 
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As before, we are now interested in describing the ground state of S^ + S^^^, which we expect 
carries total spin S = 1/2. However, now the projection onto the state with S = 1/2 cannot 
be done by the method used for the LGW theory. For g > gc, we begin with a 5 = ground 
state of S§, but now don't find that the impurity term in Eq. (11.81) introduces any net spin: 
the total spin remains at S" = to all orders in perturbation theory. Clearly, we need the 
impurity charge Q to non-perturbatively bind a 5" = 1/2 Za spinon. For g > gc, such binding 
can be addressed via a non-relativistic Schrodinger equation^, the analysis does not appear 
appropriate at the main point of interest, g = g^, where we have a conformal field theory 
(CFT) with no sharp quasiparticle excitations. Here we expect the spinon to be smeared 
over the whole system of size of L. We shall describe this spinon state by explicitly beginning 
with a 5* = 1/2 state of and then perturbatively examining the influence of iSj^p: this is 
expected to yield correlations in the true S = 1/2 ground state of + S^^^. 

Using the language of general SU(A^), let the ground states of + iSj^p be these 
transform under the fundamental representation of SU(iV). To find the matrix element of 
some operator 0{x) between states \a) and of the SU(A^) multiplet, we compute, 

(^,(0, r/2) exp (-1 /5^/2 ^-(0' Oix, 0)4(0, -T/2)\ 

(«|0(f)|/5)= hm ^ ^- (1.9) 

{za{0, T/2) exp [-1 j^^)^ A,{0, r)rfrj zl{0, -T/2))^^ 

Effectively, we start with external charge free vacuum, and then at time r = —T/2 create a 
spinon together with the Wilson line, the latter representing the effect of the external charge 
Q = 1. We wait for a long time T/2 to single out the lowest energy state with the quantum 
numbers of the operator 2^. We then measure the operator 0{x), again wait time T/2 and 
annihilate our spinon together with the external charge. The denominator in Eq. (11.91) serves 
to cancel out the matrix element for creating the spinon - external charge bound state out 
of the vacuum (no sum over a is implied in the denominator). Expressions of type (II. 9p are 
common when studying the properties of heavy-light mesons in quantum chromo dynamics. 

The time T must be much larger than the gap between states with the quantum numbers 
that we are studying. In the spin gap phase, g > g^ this gap is finite in the infinite volume 
limit. However, at the critical point the gap will be of order 1/L. So one has to choose 
T ^ L. Although unusual, this condition can always be satisfied as we work at zero 
temperature. 

To discuss higher charge impurity {Q > 1) one needs to act on the vacuum with higher 
U{1) charge composite operators of the z field. The resulting states can form higher repre- 
sentations of SU(A^) symmetry. For simplicity, we limit ourselves to Q = 1 below. 

Details of our evaluation of Eq. (II. 9p in the 1/A^ expansion appear in Section [111 Al We 
will obtain results for the scaling functions appearing in Eq. (11.31) and (II. 5p describing the 
spin distribution at the deconfined quantum critical point. 

In addition, in Section IIII Bl we compute the uniform and staggered spin distributions 
in the Neel phase of the CP^~^ model. We find that the short distance behaviour of spin 
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distributions both at the critical point and in the Neel phase is in agreement with the 
impurity scahng theory postulated in Ref. |30|. In particular, we obtain substantial additional 
evidence that the uniform and staggered spin operators flow to the same impurity spin 
operator upon approaching the impurity site. Results of the expansion for the impurity 
critical exponents of uniform and staggered magnetization are obtained. 



II. LGW CRITICALITY 

This section will study the field theory 5" + describing an impurity in an antifer- 
romagnet with an even number of S* = 1/2 spins per unit cell. As discussed in Section [U 
the 0(3) non-linear sigma model formulation in Eqs. (11. ip and (11.21) is not appropriate for 
our purposes. Instead, we shall use a 'soft-spin' approach which yields a convenient descrip- 
tion of the rotationally-invariant state of the bulk antiferromagnet for g > g^, and of its 
impurity-induced deformations. The universal results appear in an expansion in 



e= (3-d), 



(2.1) 



where d is the spatial dimensionahty. 

This dimensionality expansions allow us to compute, in principle, the universal scaling 
functions, appearing in Eqs. (II. 3p and (II. 5p . which were numerically computed recently in 



Ref. Il8|. The scaling functions clearly depend upon the geometry of the sample, and the 
nature of the finite-size boundary conditions. Such features are not easily captured in a 
dimensionality expansion. Consequently the results in this section are more a "proof of 
principle" that the scaling results apply. Direct comparison of the results below for scaling 
functions to the numerical results are not very useful. 

As discussed in Ref. llO, the e expansion is obtained by replacing the fixed length field n by 
a field (f) whose amplitude is allowed to vary freely. However, we do not have the freedom to 
relax the length constraint on the impurity site because the Berry phase term is only defined 
for a unit length field. Consequently, we retain an independent field n(r) representing the 
impurity spin, which is now linearly coupled to (p. So we consider the theory 

= J V(f>{x, T)Vn{T)6 (n^ - l) exp - 



imp 



s 



d xdr 



t [io,cf>r + s<i>^) + 1 (0^) 



2\2 



St, = ^SJ dTA[n{r)] ■ ^ - 7o5n(r) ■ 0(0, r) 



(2.2) 



Here s ~ (? is the coupling that tunes the system across the bulk quantum phase transition, 
and go and 70 are the couplings which were shown in Ref. ilO| to approach fixed point values 
in the vicinity of the quantum critical point. In the (3 — d) expansion, these fixed point 



values are small with (/o ~ 7o ~ ^- It "was argued in Ref. Il5| that this fixed point is identical 
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to that obtained from the 0(3) non-hnear sigma model theory appearing in Eqs. fll.llll.2p . 

We will be interested here in the s > Sc regime of here, where (0) = and full 
rotational symmetry is preserved in the absence of the impurity. As discussed in Section [H 
we need to project on to the state with total Sz = 1/2 in the presence of the impurity. This 
is easily done here by choosing the following parameterization for the impurity degree of 
freedom n(r) in terms of a complex scalar iP{t): 

^ + ^* TTPJ ^ - V^* 



n = (^^^^ V2 - ^^^^2 - 1 - I^I^J . (2.3) 
The advantage of the representation 02.31) is that with the gauge choice 

A(n) = -^(-n„n,,0), (2.4) 
the Berry phase takes the following form 

Furthermore, the measure term in the functional integral also has the simple form 

J Vn5 (n^ - 1) = y V'ipV^* (2.6) 

Now, an expansion of the correlators of 2^, in a functional integral over (p and ip about the 
saddle point with = and ip = 0, in powers of the couplings 70 and go, automatically 
projects onto the state with total spin projection Sz = 1/2. This is easily established by 



applying a uniform magnetic field, and verifying by the methods of Ref. llOlJlSl that the total 
magnetization is quantized by a Ward identity associated with the conservation of spin. 

We can now use the above perturbative expansion, using methods explained at length 
elsewhere^ii^, to compute the expectation values of the magnetization density {Qz{x)) and 
the Neel order parameter (0^(af)). We perform this computation on a sample with periodic 
boundary conditions and length L in each spatial dimension, i.e. a torus T'^. The main effect 
of the finite boundary conditions is that the momenta p are discrete, and each momentum 
component is quantized in integer multiples of 27r/L. The results below are easily generalized 
to other finite size geometries and boundary conditions. To leading order in e, the results 
are 



(Qzix)) = S5\x) - ^',SS\x) I ^^-L-^G{co,0) + 2^',S J ^G{co,x)G{co,-x) 



(0,(f)) = 70^^(0, f) 



(2.7) 



where e is a positive infinitesimal proportional to an applied magnetic field which selects 
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the Sz = 1/2 state. We may set e = after the frequency integrals have been performed. 
The Green's function of the (p field is 



p 



(2.8) 



where A is the spin gap of the bulk antiferromagnet in the absence of the impurity. Other 
boundary conditions will only change the form of G. It is easy to check that the spatial 
integral of (Qz) is quantized at 5". 

To leading order in e, it would appear that we can set A equal to the spin gap in the 
infinite bulk antiferromagnet, and in particular, set A = at the critical point s = Sc- 
However, we will see below that for the particular boundary conditions we are using here, 
there are infrared divergencies at A = in the expressions for the impurity-induced spin 
textures. In such a situation we have to examine the finite L corrections to the value of 
A at s = Sc, which yield a non-zero A even at the bulk quantum critical point. The value 
of A can be computed as described elsewhere^, and to leading order in e, the equation 
determining A at the quantum critical point s = Sc is 



^90 1 
6 



1 



27r u;2 + p2 + 



(2.9) 



To leading order in e, only the p = term on the right-hand-side has to be included; setting 
go equal to its fixed point valued we find for small e 



A 



11 



1/3 



(2.10) 



Note that LA is a universal number at s = Sc, which is the main result we will need below 
to establish the universality of the spin texture. 

Returning to the expressions in Eq. (12. 7p . we now want to manipulate them into the 
forms of Eq. (II. 3p and (II. 5p . However, the presence of the 6'^{x) in Eq. (12. 7p makes the x 
dependence singular. These singularities are in fact an artifact of the present perturbative 
expansion in real space, and are not expected to be present once the expansion is resummed. 
This is evident by examining the results in momentum space, where the results are a smooth 
function of momentum. In this manner we obtain after applying Eq. (12. 8p to Eq. (12. 7p 



(QM) 
{Up)) 



s 



p2 + A2 



2El 



(2.11) 
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where Ej^ = + A2. Now Eqs. (KWt can be evaluated at the fixed point value of 70, 
and to leading order in e they are seen to yield results consistent with the following scaling 
forms which can be deduced from Eqs. fll.3lll.5p 



(QM) = ^QipL) 



(2.12) 



The explicit results for the scaling functions to leading order in e are 



s 



1 - 2n'e 



y—(- 



— — \l — 27r^e (a finite number)] 

,-72 _L 1^2/^2 L n 



y 



(2.13) 



where now x and y are three dimensional momenta whose components are quantized in 
integer multiples of 2tt (except in the integral in the second equation), and £g = \/a? + L^A^. 
It is easily checked that these expressions are free of infrared and ultraviolet divergencies, 
and so yield universal results because LA is a universal number. 

From the above expression, we observe that ^gdy] — > 00) = ^(l — (e/2)ln|^), which 
we assume exponentiates to $q(|^| — >■ 00) ~ \y\~^^'^. From the short distance behavior of 
the spin texture discussed in Ref. [181, we expect that $Q(|y| — > 00) ~ \y\~'^'^'^, where rj' 
is the scaling dimension of the boundary spinr^. So we obtain the value 77' = e, which is 
consistent with earlier results^. Similarly, from the short distance behavior discussed in 



Ref. |l8|, we also have ^'n(|y| — ^ 00) 



\y\ 



-2+{e+-n-n')/2 



So with T] ~ C^(e^) and t]' 



have $„(|y| — > 00) ~ |y| ^, which is consistent with Eq. (12.131) . 



III. DECONFINED CRITICALITY 

This section describes the Neel-VBS transition in square lattice quantum antiferromag- 
nets with a single S = 1/2 per unit cell. As discussed in Section HI the response of a 
non-magnetic impurity is described by the action + S^^^ in Eqs. (11.7111. 81) for a complex 
SU(A^) spinon field Za and a non-compact U(l) gauge field A^. Here we will describe the 
1/A^ expansion of its universal critical properties. Note that in what follows we have rescaled 
the spinon field z, to remove the coupling constant g from the action (II. 7p . in favour of a 
rescaled constraint z^^Za = 1/g. This constraint is enforced with a local Lagrange multiplier 
A, so that the bulk action becomes. 



5f 



dr I d X 



(3.1) 



It is useful to define SU(A^) generalizations of the SU(2) observables introduced in Sec- 
tion [H The uniform magnetization density Q generalizes to Q"", which is the temporal 
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component of a current associated with the SU(A^) rotation symmetry, 



Q 



(3.2) 



(where = d^ — iA^ is the covariant derivative) while the Neel order n in Eq. (I1.6P becomes 
the staggered magnetization operator 



where are generators of the SU(A^) algebra. We will describe the spatial dependence 
of the expectation values of these operators for two cases: a finite system of size L at the 
critical point g = Qc iii Section 1111 At and the infinite system in the Neel phase with broken 
SU(A^) symmetry in Section fill B[ 

A. Critical point in a finite system 

We tune the system to the critical point g = Qc of the infinite volume zero temperature 
model, and then consider the system on a spatial torus of length L. We use periodic boundary 
conditions for all fields.^ As we discussed in Section HI the ground state in the absence of an 
impurity is a spin-singlet, while adding an impurity yields a ground state which transforms 
under the fundamental representation of SU(iV). This ground state has a single spinon in 
it, and we argued that the projection onto this state can be performed by Eq. (11. 9p . For an 
additional test of our projection formalism, see the appendix, where we compute the U{1) 
(electric) charge density in the presence of the impurity. 

Before we address the explicit computation of (II. 9p . we discuss scahng forms that our 
results should obey. 

1. Scaling Forms 

We are interested in computing the uniform and staggered magnetization densities. Re- 
call, that since the uniform magnetization is a zeroth component of a conserved current, it 
receives no renormalizations. Therefore, utilizing the SU(A^) symmetry, we have the general 
scaling form. 



^ In principle, on a spatial torus, we can certainly have a finite magnetic (Fij) flux, which would correspond 
to non-periodic boundary conditions. However, finite flux sectors are expected to be separated from 
vacuum by an energy gap, and hence are suppressed at T = 0. 



(3.3) 




(3.4) 
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The leading prefactor corresponds to the scahng dimension Aq = d = 2 oi the magne- 
tization density. Moreover, by conservation of total SU(A^) charge, 



rfV<l>Q(r) = -l (3.5) 

where the integral is over < ri,r2 < 1. Similarly, for the case of the staggered magnetiza- 
tion, 

iV""- fx 



{a\n\xm = A^" (^-J $„ (^-J (3.6) 

Here ?7„ is the anomalous dimension of the staggered magnetization operator n"(x), A„ = 
dim[n°'] = 1 — rjn. This exponent is related to the exponent t] in Eq. ( 11. Sp . and their 
values were computed previously^^ in the expansion for arbitrary spacetime dimension 
2 < D < 4: 

The function $q is completely universal, whereas is universal only up to an overall scale. 
In particular, does not have any property analogous to ( 13. Sp . 

Of particular interest is the behavior of the functions $Q(r), $ri(r) for f 0. We make 
a hypothesis that n°(x, r) and Q°'{x, r) flow to the same operator 5''^(r) as x approaches the 
Wilson line, 



,lim^Q'^(f,r) = -^5«(r) (3.8) 



|x|-*0 



lim n"(x,r) = " g"(r) 

|x|^0 \x\ ™P 

Calculations in the e expansion supporting this hypothesis have been given in Ref. |3d|. We 
have performed analogous calculations in the expansion again confirming the OPE (13.81) . 
Technically, this impurity OPE program consists of the following steps. First one considers 
the (multiplicative) renormalization of the operator n"(x = 0), by studying its insertion 
into the two point function of the z field (this consist of the usual bulk renormalization, 
plus an additional renormalization of the logarithmic divergences that appear as x — > 0). 
Once n"-{x = 0) operator is renormalized, one considers the insertion of Q°'{x 0) into the 
two point function of the z field. The highest divergence as \x\ — is power-like, 
modified by logarithms at higher orders in This leading divergence can be cancelled 

by a n°'{x = 0) counterterm (with a coefficient that diverges as x —>■ 0). This procedure 
gives one a way to construct order by order in the impurity operator S°'{t) (which 

is essentially a regularized n°-{x = 0,r)), and compute the anomalous dimensions A|^p, 
A[J^p as well as coefficients cq, Cn (the later are renormalization scheme dependent). As the 
computation of the OPE in the expansion essentially follows that in the e expansion 



presented in Ref. l30|, we shall not include it here. We only note that in this way, we have 
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been able to explicitly check the OPE (13.81) to order obtaining Af^^ to order 

and A|^p to order (this is lower order than the corresponding result for A]!^^ as cq/cn 

is of order 1/N). Explicit results in this expansion will appear in Section fill B[ 

Calculations of $q and given below provide additional support for the OPE (13.81) . 
Note that the exponents A|^p and AJJ^^p are not independent. Indeed, let the correlator 

{S''ir)S\0))r^-^Sa,. (3.9) 



The exponent A^ is related to the boundary spin exponent rj' used in Refs. Il0l]l8l by rj' = 2 As- 
Then, 

Ac. = Aq + Ag,p = A„ + AC„p (3.10) 
Recalling, Aq = 2, A„ = 1 - ?7„, 

A«p = ALp-l-r^„. (3.11) 

Our explicit results for the profiles $q,$„ confirm the relation (13.111) to leading (zeroth) 
order in 1/A^, see below. We have also been able to check this relation to order 1/A^ 
using the impurity OPE program summarized above: to this order, A|^p = — 1 — rjn, as 
^imp ~ 0{1/N'^). The result of our evaluation of Af^^ to 0{1/N^) will appear later in 
Eqs. ( I3:98|) . ( 13:991) . 

Note that the OPE (13.81) is sensitive only to short distance physics, and, thus, coefficients 
cq, Cn should be independent of the system size L as well as the deviation from the critical 
point (all this IR information is, however, contained in the impurity operator S""). Thus, 
the ratio, 

— = lim \x\ ™p ""p- — T— r— = lim \x\ ^1"- — — — (3.12) 
Cn \x\^o {n"-{x)) \x\^o (n"(x)) 

although non- universal, should be constant throughout the scaling regime (once the reg- 
ularization scheme is chosen). We shall check this fact below to leading order in 1/A^ by 
comparing the short distance behaviour (controlled by the OPE) of uniform and staggered 
magnetization densities at the critical point and in the Neel phase. 



2. Projection onto the Single Spinon State 

Now we return to the evaluation of the matrix elements (II. 9p . Although it is possible 
to obtain all the results presented below directly from Eq. (II. 9p it is technically somewhat 
simpler to use instead, 

(zJk, r/2)0(f , 0)zUk', -T/2))imp 
{a\0{x)\P)= lim ^ ^ ' : ^ ^ \ V (3.13) 
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Here, Zaik^r) = J (f )e and the subscript "imp" indicates that the correlator 

should be computed in a theory with the action + S^^^ which includes the impurity 
term. Effectively, we have extended the Wilson line, which in (11.91) stretched from the point 
where a spinon was created to the point where it was destroyed, to run from r = — oo to 
T = oo. In addition, we have taken our "incoming" and "outgoing" spinon to be in momenta 
k and k' states. This makes the numerator and denominator of ( 13.13^ non-gauge invariant. 
Nevertheless, we expect that this non-gauge invariance comes solely from the matrix element 
for creating the ground state of the system by acting on the vacuum with z'^ and cancels out 
between the numerator and denominator of (13.131) . 

Since the impurity term Eq. (II. 8p breaks spatial (but not temporal) translational invari- 
ance, for T — s> oo we expect to obtain the ground state irrespective of which /c, k' we started 
with. Nevertheless, it will be most convenient in our perturbative treatment to work with 
k = k' = Q. 

Since the external charge does not break SU(A^) symmetry and time translation symme- 
try, we have, 

(Za(x)2;jj(x'))imp = 5apD{x, X , T - t') (3.14) 

We let, 

D(f,f ,r) = ^,Y1 j ^^(P,p',^)e^^^e-^^^^'e- (3.15) 

We write, 

{za{y)0{x)z'^p{y'));a^p = j dvdv'D{y,v)0ai3iv,x,v')D{v',y') (3.16) 
Fourier transforming, 

p p' 

(3.17) 

where we use the notation that the three-vector x has spatial components x and temporal 
component Xr- So, 

{z^ik, r/2)o(x, Q)zl{k\ -r/2))i,,p = 

E / '^^'^-^D{k,p,u:)0^,{p,qJ,u,J^^^^ (3.18) 

^ — ' J ZTT ZTT 

As we perform the integral over tu, uj' , we pick up poles of the propagators D in the ^iyj) > 0, 
'^{oj') > planes (we expect that Oag is analytic in uj). In the limit T — oo only the 
contribution from the pole with smallest imaginary part survives. Let this pole be at = im 
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and denote by Res{k,p) the residue of D{k,p,to) at this pole. Then, 

(?Jik, T/2)0(f, 0)4{k', -T/2))i^p^5^(zRes(fc,p))(^Res(p-', k'))O^M qj, im, 2m)e*^^V'"^ 



(3.19) 



Similarly, the denominator of fl3.13p is, 

{z^{k, r/2)zi{k', -T/2))i^p ^ LHRes{k, fc')e-™^ (3.20) 

Finally, 

{a\0{xm = -l^(a|0(g-)|/5)e^«^ (3.21) 

9 



with. 



5. Large N expansion of CP^ ^ theory in finite volume 

We now compute the expression (13.221) using the large N expansion in finite volume. 
First, consider the N = oo limit. The gap equation reads, 

and to this order in A^, ml = i{X). In the infinite volume, the critical coupling g = Qc is 
obtained when the gap mo vanishes, 

' (3.24) 



g,N J (27r)3p2 

However, once we make the spatial volume finite, a non-zero mo is generated even at the 
critical point. Thus, setting g = g^, using Eq. (13.241) and poisson resumming, we obtain, 

^ I 27: J (27r)2 u;^ + ^ + J 2ti J {2Tifuj'^+f ^ ' ' 



On the left-hand side, only the n = term diverges in the UV. However, this divergence 
cancels with the divergence of the right-hand side. Thus, performing all integrals, 

y J_^-mo\H\L ^ (3.26) 
477 77, AtT 
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The solution of the Eq. (13.260 is, 

mo = 9y (3.27) 

where ^ is a constant that can be obtained by solving (13.261) numerically to be, 6 ^ 1.51196. 
Thus, at leading order the propagator, 

Do{k, k\u) = (3.28) 

and the lowest pole is at = 0, = irriQ and, 'iRes(A;, p) = ^kQ^pfl^^^- 

To develop the 1/A^ expansion, we will need to find the and A propagators. The 
dynamically generated self-energy for is to leading order. 



dQr f (2g-p)^(2g-p)^ 25^,^ . ^g) 

{{q - pY + ml){q'^ + ml) {q^ + ml) ' 



This self energy is always more singular near the critical point than the bare Maxwell term 
in iS^, and so we will work with = oo for the rest of this paper. To find the photon 
propagator, V^^{p), we also need to fix a gauge. Practically, for the calculations to follow, 
we will only need the static electromagnetic propagator P^^(p, p,- = 0) = K^T-ip^Vr = 0)^^ 
which is a gauge invariant quantity. We also note that in the infinite volume limit, 

K,M = K{q){q\,-q^q,) (3.30) 
K{q) = NAq^-^ (3.31) 

where the constant A is given by, 

^_ 1 {D-2)T{2-D/2)TiD/2-l)' 

^~J4^ fiD) ^"^-"^^^ 

Here D is the space-time dimension. In our case, D = 3 and A = j^. 
Likewise, the self-energy for A is to leading order, 

n(p) = N^Y.I ^ ,2, 2^// yT—^ (3-33) 

^ 7 2ti {q^ + ml) {{q - p)^ + ml) ^ ' 

In the infinite volume limit, 

n(p) = NBp^-"^ (3.34) 

where the constant B is given by, 

^_ 1 r(2-D/2)r(D/2-i)^ 

^ - {A7t)d/2 r{D - 2) ^^-^^^ 
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FIG. 1: The insertion of into the z propagator. 



For D = 3, 5 = \. 



4- Matrix Elements 

Now, let us compute the matrix elements of operator Q"-{x). The insertion of into the 
z propagator, to leading order in is given by diagram in Fig. [H so 

Q^^(p,g,p',w,w') = ^(^ + ^')7^:A>-'-p (3-36) 

So utilizing formula (I3.22p . with k = k' = 0, we obtain, 

= -T„V,D (3.37) 

i.e., 

{amxm = -^^T^^ (3.38) 

and the function $Q(r) = — 1, satisfies the normalization condition (13.51) . So at leading order 
in the expansion the magnetization in the presence of an impurity is spatially uniform. 
The system with the impurity simply consists of a free spinon in the zero momentum state. 
The effects of the interaction with the impurity appear only at next order in 

Similarly, for the staggered magnetization, the insertion of n"'{x) into the z propagator, 
to leading order is given by the same diagram in Fig. [1], except the cross now stands for n"". 

{P, q,^,^, ^') = Sq,p'-pTap (3.39) 

so that, 

{a\n'{q)\P) = J-6,^T^^ (3.40) 

and, 

(aK(f)|/3) = ^T,«^ (3.41) 

So the staggered magnetization at leading order in 1/A^ is also uniform, $„(r) = ^. 
Now, let's include the 1/A^ corrections. 

We will concentrate on corrections to {a\0{q)\P), for O = Q°',n"', with g 7^ (where the 
leading 0(1) term vanishes). These turn out to be much simpler to compute than corrections 
for q = 0. Moreover, for Q", we know by SU(A^) charge conservation that the N = 00 result 
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FIG. 2: corrections to z self-energy. 

f l3.37p at g = receives no further corrections. Thus, to order 

{a\0{q)\i5) iRes(0, -q)x Oa_fi{-q, <?, 0, irriQ, imo)o + iRes{q, 0)i O^fsiO, g, q, im^, imo)o 
+ iRes{0,0)oOaf3iO,q,0,imo,imo)i (3.42) 

where the subscripts 0, 1 indicate the order in to which the quantity has to be computed. 

The corrections to the z self-energy are shown in Fig. [2] (we drop A tadpole dia- 
grams). Of these only the last one couples to the impurity and, therefore, breaks translational 
invariance. So, letting, 

D{k,k',iu) = Do{k,k',iu) -^Do{k,p,uj)J:{p,^,uj)D{f,k',uj) (3.43) 

p,p' 

S(fc, k', uj) ''t^' ^2iujVrr{k - k\ 0) + 0{l/N^) (3.44) 

and 

Dik.k'^u;) -j^2zuVrr{k-k',0)^^ i^^^^ -,+0{l/N') (3.45) 

So the residue, 

1 1 

iRes(0, -q) = ^Res(g, 0) -^^V^M 0) + Oil/N^) (3.46) 

Note that at this order renormalization of the location of the pole u = im ^=°° can be 
neglected. 

The corrections to the insertion of into the z propagator are shown in Fig. [31 
Again, only the last one of these couples to the impurity and breaks translational invari- 
ance, so, 

Qlp{p, g, jT, ^, u') -2^V,^{q + p- ff, 0)T,"^ + 0{1/N') (3.47) 

Combining (13:421) . (13:^61) . flCT]) . 

{amm = - l^ko + (1 - + -^)^rriq,0)) + 0{l/N') (3.48) 



17 



FIG. 3: corrections to the insertion of Q"^ into the z propagator. 

The calculation of corrections to result fl3.40p for impurity induced staggered mag- 
netization n°'{x) proceed in the same fashion. The corrections to insertion of into the 
z propagator are given by the first two diagrams in Fig. [3] (except now the cross stands for 

insertion). None of these break translational invariance (as the last diagram in Fig. [3] is 
present only for Q"", but not for n°). Therefore, 

(c^Wm = ^ (5^,0(1 + 0{1/N)) + (1 - 6,^,o)^^V^M 0)) + 0(l/iV2) (3.49) 

Note again that in the case of {a\n°'{q)\(3) we have computed the corrections only to 
g 7^ 0. Unlike the case of uniform magnetization, here the N = 00 result for {a\n°-{q = 0)1(3) 
is expected to receive corrections. 
Thus, the scaling functions, 

1 477? 2 ^ 

$Q(f/L) = -l--^J2(^ + -^)V^r{q,0)e'''' + O{l/N') (3.50) 
$„(f/L) = l + ci + -l5^4p..(g,0)e^^^ + O(l/iV2) (3.51) 

where ci is an x- independent constant of order 1/A^ (ci should be also independent of A; we 
have not verified this fact as we did not compute the 1/A^ corrections to {a\n"-{q = 0)|/3)). 
We may write, 

$Q(f) = -(1 + 1/Q(f)) + 0{1/N') (3.52) 
*n(r) = ^(1 + 2c,0 + l/„(r)) + 0{1/N') (3.53) 

We have evaluated the functions /g, numerically and plotted them along the diagonal of 
our spatial torus in Fig. HI 

Now, we would like to find the g — > cxd, x — asymptotes of fl3.50p . fl3.5ip . For this 
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Staggered Magnetization Profile 
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FIG. 4: Uniform (a) and Staggered (b) magnetization distribution functions fqir), fn{r) plotted 
along the torus diagonal. 



purpose, we may replace the finite box propagator 2^rr('?) by the infinite box propagator, 



1 1 

naJ^ 



Writing, $Q,n(x/L) = ^^ff^Q,niq)e 



iqx 



g— ►oo 



q—*oo 



1 L 



2 1 

NAL\^^ 



+ 0{l/N^' 



Fourier transforming. 



2TTdNA Ifl 



+ 0{l/N^ 



^n{^ ^ + C2 + 0(l/iV2) 



where C2 is a constant of order 1/A^. 
Thus, we conclude that. 



AP.p = -l+0(l/iV) Al^ = 0(l/iV^ 



imp 



(3.54) 



(3.55) 
(3.56) 



(3.57) 
(3.58) 



(3.59) 



which is consistent with the relation between impurity exponents (13.111) . Note that the 
present calculation shows that is zero to order 1/A^. We shall verify this fact in a 

different way in Section UlI Bl and compute to order 
Moreover, the ratio, 



cq 



TcNA 



+ 0{l/N^) 



(3.60) 
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is independent of regularization at this order in N. 



B. Neel Phase 

In this section, we compute the uniform and staggered magnetization in the presence of 
an impurity of charge Q in the symmetry broken phase, g < Qc- We work in infinite volume. 
We develop the 1/N expansion around the symmetry broken vacuum, 

(zi) = -^v (3.61) 

Note that in general v is not a gauge invariant quantity. However, this fact does not manifest 
itself at the order at which we are working. To leading order in A^, 

-v' = --- (3.62) 

Note that ~ O(A^). Moreover, we take Q ~ 0(1) in A^. 

We now must quantize our theory around the symmetry broken state. We write, 

zi = —j={h + V + Za = tTq,, a = 2 .. N (3.63) 

v2 

We work in the so-called gauge, in which the mixing between the goldstone and the 
photon is absent, at the expense of introducing a ghost field c. In what follows, we 
have eliminated the mixing only to leading order in 1/N. This is achieved by using the 
gauge-fixing condition, 

d^A^ = ^vK~'<f) + w (3.64) 

where the action for the auxiliary field w, which appears in the Fadeev-Poppov formalism, 
is 

5^ = y dxdyw{x)K{x - y)w{y) (3.65) 

Here, K{x — y) is the photon polarization function given by Eq. (I3.3ip . Similarly, in what 
follows H(a; — y) is the A self energy given by Eq. (13.331) . 
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At the end of the day, the action one obtains is, 
= \ j dxdyA^{x) (^^y{x - y) - ^d^d^K{x - y) + d^^v'^^ 

+ dxdy4){x) {-d^6{x -y) +^v^K'\x -y))(f){y) + ^ j dxdy\{x)Yi{x - y)\{y) 
+ j dxdy(^c (^—d'^6{x — y) + ^v'^K^^{x — y)) c{y) + ^vc{x)K^^{x — y)h{y)c{y)^ 
+ jdx (^\D^n\^ + ^{d^hf + ivXh + {(Pd^h - d^(Ph)A^ + {vh + + \cl>^)Al^ 
+ j dx + + h"^)^ (3.66) 

As usual, we avoid double counting by dropping any diagrams, which are already included 
in the dynamically generated N = oo self-energies for A^, A etc. The propogators for our 
fields are shown in Fig. [5l Note that in the Neel phase, we get mixing between the A and h 
fields. 

A '■wwv^w^^wwwv. A T) ir>\ — 1 ( s -^ (g-i)-^(p)PMP.^ 

^A* ^^^^^^^^^^^^^^^^ ^f^vyP) — p2K(p)+v'^ p^K(p)+(v'^ 

V\AAAAAAAAAAAAAAAA D^{p) = ^^^^J^-^^^-^ 



7r„ 



h Dh{p) 



n(p) 



p^Il{p)+v'^ 



2 

\ \ Dx{p) = p2n(p)+t,2 

h A Dhx{p) = p2n(p7+i«2 

FIG. 5: Propagators in the Neel phase. 

Now, having set up the perturbation theory, we wish to compute, {Q"-{x)), {n"-{x)). 
Utilizing the pattern of spontaneous symmetry breaking, U{N) —>■ U {N — 1) (here we look 
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FIG. 6: Leading contribution to uniform magnetization in the symmetry broken phase. 



only at global symmetry), one can show that, 

(n'^)=T{\(n°) (3.67) 

where n° = z'^T^z and T° is any generator of SU(A^) with T^^^ = 1. Similarly for Q"^. For 
definiteness, we may choose T{\ = 1, T^^ = Tj^-^ = 0, = —j^Sap, a, P = 2 .. N. 
Let's start with computing the uniform magnetization. 

Q° = ^^jl —Jr (3.68) 

where 

jI = zlDrZ, - {D,z,yz, (3.69) 

and jr is the U{1) charge density discussed in the appendix, see Eq. ( lAll) . By equation of 
motion flA4l) . 

(jV(f)) = -Jr\x) = -Q6'{x) (3.70) 
So, it remains to compute (j^(s)). Expanding jl in terms of (j), h and A^, 

jl = -iv^Ar + iv{dr<p - 2Arh) + i{hdr<p - (pdrh - Ar{h^ + 02)) (3.71) 

In the expansion the leading contribution to {jl) is of 0(1) and comes from the first 
term on the r.h.s of (13.711) . see Fig. [61 

= -Qv^VrAp, 0) = (3.72) 

Thus, (Q°(p)) = {ilip}) +0{1/N). Fourier transforming, 

mS)) ^ ^^AHjm-vSm)^ OiHN) , (3.73) 



2tiNA \x\ 4A^M2 V V / V , , 

where Hq is the Struve function and Yq is the Bessel function. Taking the short and long 
distance asymptotes. 
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FIG. 7: Leading x-dependent contribution to staggered magnetization in the symmetry broken 
phase. 

The long distance decay is a consequence of the Goldstone physics of the spin waves, and the 
decay is expected to be exact. At short distances, we have the physics of the critical 
point, and the exponent will have corrections at higher order. From the present result we 
can conclude that the impurity exponent 

A?mp = -1 + 0(VA^), (3.76) 

which is consistent with the result obtained at the critical point (13.591) . 
Now, let's discuss the staggered magnetization, 

n° = -J^z\zi - ^^^—z'^z (3.77) 



By equations of motion, 

thus, 

and 



A^-l ' N-l 



z^z = - (3.78) 



n° = ^^z]zi - ^— (3.79) 

zlz^ = + vh + hh^ + (f)^) (3.80) 



2 2 
Thus, at leading order, {zlzi{x)) = and 

{n\x))=^-v' + 0{l) (3.81) 

Moreover, the x-dependent corrections to {n^{x)) come only at 0{1/N), with diagrams of 
Fig. [7] (the part of ?7.° which contributes at this order, denoted by x, is vh). 
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a) 

FIG. 8: Leading contribution to the three point vertex of ^4^, Ai, and A fields, T^^'^(q,p,q — p). 

We will discuss the diagrams in Fig. [7] shortly. For now, we can conclude that. 

Alp = 0{1/N') (3.82) 



in agreement with the result 03.591) obtained at the critical point. Moreover, we can now 
compute the ratio, 

^ +0{1/N^) (3.83) 



Cn rcNA 

which exactly agrees with the result obtained at the critical point (13.601) for Q = 1. Notice, 
that this is a highly nontrivial check of the OPE (13.81) as {Q""), {n"-) depend on v in the 
Neel phase and on L at the critical point. Nevertheless, all the dependence on the IR scale 
cancels out in the ratio cq/c„, which is constant throughout the scaling regime. 
Coming back to the diagrams in Fig. U\ 

{n\q)) Q'v'D.iq, 0) J -^^^ {^ + \^~'^^^ ^)^^^^^^ ^^P' 0' 0))^^r.(p, Oprr{q-p, 0) 

(3.84) 

We keep the space-time dimension D arbitrary in what follows, as we wish to compare our 
result for obtained in the expansion, with the result obtained using e expansion^. 
Here, T^'^{q,p,q — p) is the lowest order contribution to the A^, Ay, A vertex, given by the 
sum of the loops in Fig. [8l The diagram in Fig. [8] a) is given by, 

rr(g,p,g-p) = 2z5^,n(g) (3.85) 

Thus, diagrams in Fig. 7 a) and b) cancel (by the way, these diagrams are individually UV 
divergent for D <?>). So, calling the diagram in Fig. [H]b), T2^{q,p,q —p), 



{n'm'=Q'v'D,{q,0) 



jD-I 



i2n 



• 0, P, 0, q- p, 0)Vrr{p, 0)Vrr{q- P, 0) 

(3.86) 
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Evaluating T'^" , 

d^^l {2l-p)^{2l-p-q), 



rr(97P,g-p) = -2iN j 



(27r)^ P{l-pY{l-qY 



(4 - D)i2x^q + (2x2 - l)p);.((2xi - l)q + (2x2 - l)p)u 



(3.87) 

where, 

= xi(l — + ^2(1 — X2)p'^ — 2X1X2P ■ q (3.88) 
We are interested only in FY, with p^ = q^ = 0. Thus, 

r7(g, 0,p, 0, q- p, 0) = - ^'^|]^^^/f y c/xic/x2C?X35(1-Xi-X2-X3)(A2)^/2-2 

(3.89) 

For IpI ^ r2'^(g, 0,p, 0, g — p, 0) ~ so for p ^ cxd the integrand in Eq. fl3.86p 

behaves as \p\~^ and the integral is UV convergent. 

We now attempt to understand the behaviour of fl3.86p for q ^ 00, from which we should 
be able to extract the impurity anomalous dimension Aj^p. For this purpose, we may set 
t> = in the propagators Dh{q,0), VrT{p,0), Vrriq — P,0) (this does not introduce any IR 
divergences) . 

<"°®> mI'T"/ (S^^rrteO.P.O.f-p.Oj^i^-lj^ (3.90) 

Let us first discuss the limit D = 4 — e, e —>■ 0. In this regime, to leading order in e, 

TYiq, 0, p, 0, q- p, 0) = -2ziV^r(2 - D/2) (3.91) 



and, 

Fourier transforming. 



^ = ^2 n?3 (3-92) 



{n\x)) + C3 - ^^^^^vHogiv\x\) + C4 + Oil/N') (3.93) 

2 iV^ 

where C3, C4 do not depend on x and are of order 1 and respectively. Thus, to leading 
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order in e, 

iV2 



Alp = ^ (3.94) 



in agreement with the calculations of Ref. l30|, where the impurity exponents were obtained 
by performing the impurity operator renormalizations as summarized in Section flll A 11 ^ 

For arbitrary AJJ^p is difficult to calculate analytically, as F"^"^ is no longer a constant. 
However, combining Eqs. (I3.86p . (13.891) and introducing a new set of Feynman parameters, 

where the numerical constant f{D) is given by, 

1 



A25(47r)^-iF(D/2- 1)2 



/(^) = A-?-n/A ^n-^T^/^,lr, 7T? / ^Xi / 6x2 / dyi / dy. 



JQ Jo Jo 



n 



(0:2(1 - X2)^?/2(l - 2/2) + a;i?/i((l - xi){\ - X2) - 2^23:2(1 - X2) - yiXiX2))~^ 

(3.96) 



Consequently, 

-S» + e3 - (4,),.-.var((J-l)/2) ^'^'l^''^ log(„^/«--|f |) + c. + 0(l/iV^) 

(3.97) 

and 



Evaluating f{D) numerically for D = 3, 



Q 



2 



ALp^-25.9^ + 0(l/iV=^) (3.99) 

We note that we have separately verified the result (13.981) by performing the impurity OPE 
program as summarized in Section IIII A II 



^ Note that in the e expansion of Ref. I,30| only the analogue of the diagram in Fig. [7] a) appears, while 
the diagrams in Figs. [7]b), c) do not appear at leading order in e, as they are higher order in coupling 
constant. Nevertheless in the 1/iV expansion, we saw that the answer comes entirely from the diagram in 
Fig. [7]c), with diagrams in Fig. [7] a) and Fig. [7|b) canceling for all D. The reason is the following: in the 
expansion all diagrams in Fig. [7] are individually of same order in e. Moreover, to leading order in e, 
the diagrams b) and c) cancel, so a) = - b) c). In the e expansion, this fact is foreseen in advance: the 
1 /e pole must cancel between diagrams b) and c) (the 4-point diagram with two photons and two scalars 
is not divergent). Thus, we can obtain the answer to leading order in e cither from a) alone or from c) 
alone. 
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IV. CONCLUSIONS 



A recent numerical studj*^ examined the spin distribution in the vicinity of a non- 
magnetic impurity in a double- layer, 5 = 1/2 square lattice antiferromagnet at its quantum 
critical point. The ground state of the system has total spin S = 1/2, and the spin dis- 
tribution of this S" = 1/2 was found to be extended across the entire system. Universal 
scaling forms (Eqs. (11.31) and (11.51) ) for the uniform and staggered spin distributions were 
postulatec^^, and found to be in excellent agreement with the numerical results. 

This paper has presented the field-theoretic foundation of the above results. Using the 
soft-spin 0(3) LGW field theory in Eq. (12. 2p . we found that the universal scaling forms in 
Eqs. (11.31) and (II. 5p were indeed obeyed in an expansion in (3 — d) (where d is the spatial 
dimensionality), and explicit results for the universal scaling functions appear in Eq. (12.131) . 

Next, we examined a similar non- magnetic impurity in S" = 1/2 antiferromagnets which 
have a single 5 = 1/2 spin per unit cell. Such antiferromagnets can display a deconfined 
quantum phase transition^iiS^ between Neel and valence bond solid (VBS) states. An explicit 



example of such a transition was found recently in Ref. 1231. We expect that such studies 
will be extended to include non-magnetic impurities in the future, and so have provided our 
theoretical predictions here. The field theory for this situation in + <Sijnp in Eqs. (11.7111.81) . 
It describes the dynamics of a SU(A^) spinor field, Za (the spinon), and we obtained its 
critical properties in a 1/A^ expansion. Projecting onto the total spin 5 = 1/2 sector 
of this theory (which contains the ground state in the presence of the impurity) was not 
straightforward here, and we achieved this by the relation Eq. (11.91) . Our results obey scaling 
forms which appear in Section [III A 1[ The scaling functions are in Eqs. (13.501) and (13.511) . 
and are plotted in Fig. |H The boundary spin exponent for the deconfined critical point 
appears in Eqs. (13.981) . (I3.99p . We also obtained substantial evidence for the structure of the 
operator product expansion near the impurity, and the fact that the staggered and uniform 
magnetizations flow to the same impurity spin operator. 

Our study of the deconfined case has so far only examined the uniform and staggered spin 
magnetizations near the non-magnetic impurity. We have not provided here a description of 
the structure of the VBS order near the impurity. This will be presented in a forthcoming 
paper. 
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APPENDIX A: U{1) CHARGE DENSITY 

Throughout the paper we have concentrated on computing matrix elements of uniform 
and staggered magnetization Q°'{x), n"(x). However, for the deconfined critical point, it is 
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also interesting to compute the charge density associated with the U{1) local symmetry of 



the CP^ ^ model. This charge density is the zeroth component of the current, 



Zj'Z 



(Al) 



As we shall see this computation serves as an additional test of our procedure for projecting 
onto the single spinon state. 

Consider the CP^^^ model coupled to an external current, 



(A2) 



As in the rest of the paper, we set = cxo, so that the gauge field has no bare kinetic term. 
Then, by equations of motion. 







5S_ 



J„ 



(A3) 



Jm = -JT (A4) 

Thus, the dynamical current completely screens (locally!) the external current. Eq. (]A4p is 
an operator identity, and should, in particular, hold in the ground state of the system with 
a single impurity. Let's check this statement in the expansion. 
We start from Eq. (11. 9p . with 0{x) = jo{x). We write. 



{a\jo{x)\p) = p{x)5ap 



(A5) 



with p[x) = -J2 YlqPi^^^'^^ ■ "^^^ Wilson line term in Eq. (11.91) can be incorporated into the 
action as coupling to an external current, J^^{x,t) = 5^q5^{x)6{T /2 — t)0{t + T/2). At 
leading order in 1/iV the numerator of Eq. (11.91) is given by diagrams shown in Fig. [9l while 
the denominator is given by the bare propagator D{x = 0, T). 




+ 




+ 



FIG. 9: Diagrams contributing to U{1) charge induced 

Thus, we can distinguish two contributions: the disconnected one, p^{q), coming from 
the first line in Fig. [9l and the connected one, p^{q), coming from the second line. We note 
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that the scalar loops contributing to p^{q) are precisely the same as those contributing to 
the self-energy of field, thus, 

p\q) = - j ^KoM qr)V.x{q, qr)JT\q, Qr) (A6) 



Now, 



Thus, 



K,Ml^,xiq) = - ^ (A7) 



Noting, J-*(g,g,)^^°"27r5(g.), 



/(g-) = -(1-5,^,0) (A9) 



For the q = part, the order of the limits g — > 0, — > is very important. In our finite 
system the q = mode is isolated, and, moreover, in our present treatment the Wilson line is 
of finite length, so we must take the q = limit first and then g^ —>■ 0. Hence, p^(g = 0) = 0. 
This is not surprising. In perturbation theory, we start with the vacuum which has charge 
0. Unless we manually project the system into a finite charge subspace (as we do in our 
treatment by acting on the vacuum with z, z'^ operators), we will never be able to see global 
screening of charge. Since the diagrams contributing to p^{q) are disconnected from the 
external z line, p^(g = 0) = 0. 

Now, the connected contribution, simply gives the charge density of one spinon in the 
/c = state, 

p\q) = -%,o (AlO) 
Putting the two contributions together, 

p{q) = 1 (All) 
p(f) = 5\x) (A12) 



as expected by equations of motion (lA4p . 

Thus, we have been able to check exact screening of external charge, which follows from 
equation of motion ( lA4p . to leading order in 1/A^. We see that local and global parts of the 
screening charge come from very different Feynman diagrams. 
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